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1. Introduction 

The search for Majorana quasiparticles in condensed matter systems is a subject of 
great fundamental and practical interest. Majorana quasiparticles have been recently 
proposed to exist in various systems which typically combine strong spin-orbit coupling 
(SOC), the proximity to an s-wave superconductor and magnetic elements [H El |3l H] . 

In particular, Majorana bound states (MBS) have been proposed to emerge in 
vortex cores on the surface of strong three-dimensional (3D) topological insulators 
(TIs) [5] or at the edges of two-dimensional TIs when brought to the proximity of an s- 
wave superconductor and a magnet P,[7]. Characteristic features of Majorana fermions 
have also been predicted in spin-orbit coupled two-dimensional semiconductors [H |9] 
and one- dimensional spin-orbit coupled semiconductor nanowires [10], and very recent 
experiments show signatures of MBS which are currently being scrutinized [HI HSj [13] . 
In these semiconductor systems, time-reversal symmetry is broken by a Zeeman splitting 
which, in addition to the proximity of an s-wave superconductor, is necessary to induce 
a TSC. In such systems, a large SOC is desireable as the vulnerability of the energy 
gap to disorder typically depends on the size of SOC [Ml [15]. If the Zeeman splitting 
is induced by external magnetic fields, low field strengths and therefore small orbital 
effects are beneficial for the superconducting correlations, which favors semiconducting 
materials with large effective g-factors. 

In this work we show that a doped HgTe-QW can implement a TSC, when s- 
wave superconductivity and a Zeeman-field is induced in the QW. In order to describe 
the HgTe-QW, we use the four-band BHZ model [T6l [T71 [18] valid near the transition 
between the inverted and non-inverted regime. In the inverted regime and in the absence 
of superconductivity and a Zeeman gap, the well-known helical edge states appear in 
the spectral bulk gap [TH [191 ISO]- Importantly, HgTe-QWs are intrinsically subjected 
to inversion symmetry breaking SOC due to bulk inversion asymmetry (BIA) [181 l2T] 
which breaks the Kramer's spin- degeneracy. We show that BIA-SOC alone is sufficient 
to induce a TSC when superconductivity and a Zeeman gap are added. However, tunable 
(extrinsic) inversion symmetry breaking terms provided by Rashba-SOC can be added 
as well. 

It is known that HgTe-QWs are characterized by a very large Rashba-SOC [22], 
which can be one order of magnitude larger compared to other III-V compound 
semiconductors [221 ES]. In [22], Rashba-splittings in inverted HgTe-QWs as large as 
30 meV were reported. In addition, the effective Zeeman gap in a perpendicular magnetic 
field was recently calculated and measured in zero gap samples to be consistent with an 
effective g-factor g* ~ 55.5 [21]. This shows that HgTe-based QWs could be a feasible 
candidate for a proximity induced TSC. 

In addition to a potentially feasible implementation of already existing ideas, we 
show that the combination of BIA and Rashba SOC provides an interesting handle 
on the topological properties of the TSC. In particular, we show that (consistent with 
topological invariants) the chirality of Majorana edge states can be reversed by changing 
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the relative strength of the two SOC terms. 

In contrast to the proposal of using ID helical edge states as a platform for 
MBS |5j, we would like to point out in this work the possibility of switching between 
chiral (Majorana) and the well-known helical edge states on the same one-dimensional 
boundary of a HgTe-QW. 

The paper is organized as follows. In section [2] we introduce the setup and the 
relevant model Hamiltonian. Subsequently, we describe the bulk properties of the 
system, emphasizing the topologically non-trivial and trivial phases of the system 
in section |3l In section HJ we solve the chiral Majorana edge states at a domain wall 
between topologically non-trivial and trivial phases and calculate the group velocity 
along the edge. Finally, we compare the chiral Majorana edge states that emerge when 
the system is in a TSC phase to the helical edge states that emerge in a TI phase 
in section |5l 

2. The Model 

We consider a doped HgTe-QW in the xy-plane in contact to an s-wave superconductor, 
as sketched in figure [H There are two inversioij|l-asymmetry induced SOC effects: 
Firstly, BIA-induced SOC is intrinsically present as a consequence of the zinc blende 
structure of the underlying crystal. Secondly, so-called Rashba-SOC arises when an 
external electric field is applied in z-direction. In addition, a Zeeman field is present 
in 2;-direction which can be induced by the proximity to a ferromagnetic insulator (via 
spin-exchange interaction), or intrinsically by using Mn-doped (Hgi-yMnyTe) QWs and 
a small external magnetic field polarising the Mn-spins [25]. Although a magnetic field 
induces also an orbital effect not taken into account here, we note that the vector 
potential of a perpendicular magnetic field induces an additional Zeeman splitting 
between the Kramer's spins with an effective giant g-factor of g* ~ 55.5 |23] so that 
only small external magnetic fields below IT would be needed. 
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Figure 1. HgTe-QW in proximity to an s-wave bulk superconductor (SC). 
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The resulting QW system is described in terms of the BHZ model 
complemented with extra terms accounting for BIA [T8], Rashba [17], Zeeman 
field fl8[ [25] . and the proximity to an s-wave superconductor [26]. The Hamiltonian 
can be written in the basis order i^Ei+y i^Hi+, i'Ei-,^Hi-,i^Ei-,i^m-, -i'Ei+, -i^m+ in 
terms of products of the Pauli matrices s, q_, r and the identity matrix, where s acts in 
pseudospin {El, ifl)-space, a acts in Kramer's spin spac^ and r acts in electron-hole 
space. The extended BHZ-model Hamiltonian reads 

H = (^e(P) + M{k^)sz + Ak^s^a;, - AkySy 

+h (cos 26*0"^ + sin 26ax) SyTz 



+ 1) 

2 

+ {B+ + B_Sz) + (A+ + A_s,) r, 



+-Ro^t; I K(yy - kya^ ) Tz 



(1) 



with 

e(P) = C- D{kl + kl\ M{k') = M- B{kl + fcj) (2) 
5± = {Be ± Bh) /2, A± = (As ± A^) /2. (3) 

Here, k^ and ky denote the differential operators —idx and —idy. 

The material parameters A, B, D, M depend on the geometry of the QW. In 
particular, the mass parameter M is tunable by the QW thickness. The BHZ-model 
was successfully used to explain experimental transport data for mass-parameters in the 
range —24 meV < M < 10 meV [27]. The parameter C defines the doping of the QW 
(tunable by top and/or bottom gates). 

The two additional SOC terms have coefficients denoted by h for the BIA term and 
Ro for the linear Rashba ternJJJ]. As the BIA term is not invariant under rotation in the 
xy plane, the angle 9 is introduced to describe the angle between the x-axis and the 
[100] crystal direction. 

The pairing terms A e and describe the effect of the s-wave superconductor on 
the El and HI bands. We neglect a possible off-diagonal proximity coupling between 
El and HI bands and without loss of generality assume that the pairing terms A^ and 
Ah in ([1]) are real |26j . 

Typical values of the above mentioned parameters are presented in table [H As 
Zeeman energy of the El and HI bands, we use Be = 1.3 meV and Bh = — 0.07meV, 
respectively. In addition, £ = 3mV/nm is assumed as an effective electric field 
invoking the Rashba-SOC. Unless stated otherwise, we also use a doping parameter 
of C = 7.8 meV and ^ = 0, as well as A^ = 0.5 meV and Ah = 0.5 meV for the 
superconducting proximity term^. 

§ The quantum number that transforms under time reversal The electron spin is not a good quantum 
number in the presence of spin-orbit interactions. 
II Rashba terms of higher order in k are neglected. 

% In reality, A^; and Ah will not be exactly equal, however, in our parameter regime C w — M (see 
below), the dominant contribution will come from A^;. 
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Table 1. Typical values for the parameters of the inverted HgTe-QW as taken 
from [UllIT]. 

3. Bulk Physics 

BIA-SOC is sufficient 

Firstly, we show that BIA-SOC is sufficient in order to overcome the 'fermion doubling 
problem' |3] by freezing out half the degrees of freedom (unless k = 0), which otherwise 
would make the spectrum doubly degenerate. Under periodic boundary conditions, kr^ 
and ky are good quantum numbers and the energy dispersion for the standard BHZ 
Hamiltonian supplemented with the BIA term can be calculated to 

E{k^, ky) = e{P) ± ^{A\k\±hy + M{k^y. (4) 

Clearly, the spectrum is non-degenerate for \k\ ^ 0. By additionally introducing the 
Zeeman terms, the time reversal symmetry is broken and thus the Kramer's degeneracy 
at \k\ = is lifted. 

The combined effect of BIA and Zeeman terms on the bulk energy dispersion is visualized 
in figure [2] (left). In the following, we choose the Fermi energy Ep tohe zero so that the 
Fermi level is situated within the gap between the ii^l-like bands opened by the Zeeman 
splitting. The doping parameter is tuned in such a way that C ~ —M. As a result, 
only a single pair of Fermi points exists and thus the system has effectively spinless 
dynamics. 
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Figure 2. Left: Bulk energy dispersion of a HgTe-QW with doping parameter 
C = 7.8 meV in the presence of Zeeman terms and BIA-SOC; Rashba-SOC is absent 
and only electron states are shown. Right: The same parameters are used but the 
s-wave superconducting proximity terms are included. The colors indicate the charge 
expectation value that ranges between e(red) and -e(blue). 



When, in addition, the proximity of the s-wave superconductor is taken into 
account, electrons with energy around the Fermi level and momenta k and —k couple. 
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opening a gap around the Fermi level as shown in figure [2] (right). The special conditions 
under which this gap closes will be discussed in the following. The Fermi points at which 
this happens are the key to the emergence of the chiral Majorana edge states [28] . 

Symmetry properties 

Before calculating the Fermi points we will briefiy review the symmetry properties of the 
two £'l-like bands closest to the Fermi level. Their bulk energy spectrum is presented 
in figure [3] in the presence of different SOC terms: BIA or Rashba or the cooperation of 
BIA and Rashba. Interestingly, although the spectrum is rotationally symmetric when 
either Rashba or BIA are present, this rotational symmetry is broken when Rashba and 
BIA are simultaneously present. 

This can be understood as follows [29]: While the effective magnetic fields due to 
Rashba-SOC are oriented clockwise perpendicular to the wave vector k, the effective 
magnetic fields due to BIA spin splitting have a tetrahedral symmetry, as shown 
in figure [3l When both BIA and Rashba are present, the effective magnetic fields 
add vectorially, decreasing or enhancing each other. The energy of the lowest-energy 
band is very sensitive to the effective /c-dependent in-plane magnetic field. Therefore, 
it exhibits a tetrahedral symmetry when BIA and Rashba interfere with each other. 

Fermi points 

In the following, we will identify the Fermi points of the spectrum. To simplify analytical 
calculations, we project the full Hamiltonian ([T]) onto an effective Hamiltonian. When M 
defines the dominant energy scale and the doping is chosen as C ~ — M, four of the eight 
bands become irrelevant. Using quasi-degenerate perturbation theory [29j ("Lowdin 
partitioning" ) , the eight-band Hamiltonian ([T]) can be reduced to an effective four-band 
Hamiltonian acting on the El-\ike electron and hole states only (see Appendix), 

H = (x{P)tz + Beo-^ + Aet^ + RohayT^ - R^kya^r^ 

(p) 

+hkx (cos 29ax — sin 29ay) — hky (cos 26(Ty + sin 29ax) t^. 

where the new parameters /i(fc^) = e(A;^) + Miji^) + (A^/c^ + /i^)/(2M) and h = Ah/M 
have been introduced. 

In the absence of BIA-SOC, the reduced form ([5]) is analogous to the effective 
Hamiltonians of different semiconductor systems, which have already been shown to host 
chiral Majorana edge states [HI EO, [ID]- Note that the Rashba energy Ur = 2m* [H] 
is kept quite small here by intention in order to study the combined effect of BIA and 
Rashba-SOC: for the values in table 1 it amounts to Ur = Rq/{^ — D — B) ^ 0.3 meV. 
On the other hand, the BIA-SOC leads to an effective Dresselhaus-SOC in the reduced 
model ([5D, which amounts to Ud = (Ah/Mf/i^ - D - B) ^ O.TmeV. In [3l], electric 
fields up to £^ ~ lOOmV/nm were estimated to be possible, which could enhance the 
Rashba energy in principle by a factor of 30, without affecting the effective Dresselhaus 
energy. 
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The effective Hamiltonian 
crossings at 



yields (for 6 = 7i/4 and A^; 7^ 0) zero energy level 



{Ro - hfkl + (i?o + h)^kl = - [ \Ae\ ± 



jj,{k 



2\2 



(6) 



This implies that for non-vanishing SOC, Fermi points exist (for B^^ — ^{k'^Y > 0) ^i^t 



k^ = {},ky = {) and \Ae\ = ^Bl - /i^ (7) 
= 0,kyy^ 0, Ro = -h and ([6]) (8) 
k^ y^O,ky = 0,Ro = h and §^ (9) 

where /iq = yu(fc^ = 0) = C + M + h'^/{2M). While the Fermi point d?]) exists also 
in other semiconductor systems that have been proposed to host chiral Majorana edge 
states [HI EQI [To] , the Fermi points ([H]) and ([9]) are characteristic of a system with two 
different SOC terms. 
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Figure 3. Bulk energy dispersion when a) only the Rashba term is included b) only 
the BIA term is included c) both terms are included. Top: energy dispersion E{kx) 
where fcj, = (dashed line) compared to energy dispersion E„iin{kx) where ky is such 
that the energy becomes minimal (solid line). Bottom: energy dispersion E{kx, ky) for 
the minimal positive energy band. Vectors illustrate in-plane spin expectation values 
projected onto the El-band {<7x{l + Sz)/2) and {ay{l + Sz)/^). 



Bulk Topological Invariant 

In this section, we identify the topological phases of the non-trivial momentum space 
topology of Hamiltonian ([1]). The fundamental Hamiltonian ([1]) belongs to the Cartan- 
Altland-Zirnbauer symmetry class D, which in two spatial dimensions is characterized 
by a topological invariant G Z [32j. The first Chern number which is commonly used 
as the topological invariant can be obtained for a non-interacting Hamiltonian with 
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non-degenerate spectrum as the integral of the Berry curvature over the first Brillouin 
zone 



where the Berry curvature J-" = dk^Aky —dkyAk^ is obtained from the the Bloch functions 
via = —'iJ2j{'^jW\9ka\uj{k)), o. = x,y, the index j running over the occupied 
states. 

For non-vanishing spin-orbit interactions 0, the first Chern number of the 
Hamiltonian ([1]) is given by 



where 9(x) is the Heaviside step function. This was verified by computing the 
topological invariant from a discretized version [33] of ( ITOl) for a regularized versioE0 
of Hamiltonian ([T]) numerically using MATLAB at several illustrative values of the 
parameters on each side of the topological transitions. Equation ffTTj) is valid when the 
BIA parameter h is small compared to the mass parameter M. 

While the general form of (ITTl) is analogous to known topological invariants of 
similar systems such as spin-obit coupled quantum wires pLOJ or 2D spin-orbit coupled 
semiconductors |[8|, this constitutes to our knowledge the first such calculation for a 
HgTe-QW in proximity to a superconductor. Importantly, this is also the first time that 
the BIA term h is included. When both Rashba and BIA are present, the sign of the 
topological invariant Ci can be reversed due to the interplay between both SOCs. This 
is similar to the combined effect of Rashba and Dresselhaus SOCs that was described 
for non-centrosymmetric superconductors |34j . 

In the absence of BIA, the relevant parameters for a topological transition to occur 
are the proximity coupling A^;, the Zeeman energy Be and the sum of mass gap and 
doping, /io = C + M. When BIA is taken into account, the topological transition is 
shifted towards higher doping as described by the term fj.o = C + M + ^ (in the limit 
that M is the dominant energy scale). The topological transition is presented as a 
function of Zeeman energy Be and doping C in figure HI 

4. Majorana Edge States 

In this section, we discuss the appearance of Majorana edge states in this system. 
With an edge we mean that there is a domain wall in the shape of a straight edge where 
one side is topologically nontrivial (Ci = ±1) and the other side is topologically trivial 
(Ci = 0). We distinguish between two cases: either (i) the topologically nontrivial region 
is vacuum or (ii) a topologically trivial QW. In the following, these cases will be referred 
as (i) hard-wall boundary conditions and (ii) soft-wall boundary conditions. Soft-wall 

+ i. e. BIA and Rashba should not both be zero 

* i.e. replacing + ^ {2 — cos{akx) — cos{aky)), k^ ^ ^ sm{akx), ky ^ ^ sm{aky) where a is a 
lattice constant that is chosen in such a way that no artificial level crossings occur. 
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Figure 4. Topological phase diagram in the parameter region Zeeman energy Be and 
doping C. The topological transition described approximatively by (blue solid line 
when BIA is present, red dashed line when BIA is absent) is compared to numerical 
calculations based on ([T]) (blue triangle when BIA is present, red square when BIA 
is absent). When BIA asymmetry is taken into account, the topological transition is 
shifted towards higher doping. The other parameters are given in table [T] 

boundary conditions (two different topological phases in a QW) can be achieved by 
varying the doping parameter C accordingly across the sample, as shown in figure [5] (b). 

Without loss of generality, we assume that the domain wall lies parallel to the x- 
axis and that the region y < is topologically non-trivial while y > is topologically 
trivial. Due to translational invariance along the x-axis we can consider solutions that 
propagate along the edge with definite k^. 

The general strategy when solving such a boundary problem is as follows: In a first 
step, we use the Ansatz 

i:{x,y,E) = e"'^^e''yyij{K,ky), (12) 

to calculate the wave functions on both sides of the boundary. Hereby, ky, which is not 
a good quantum number, has to be determined as a function of k^ and the energy E 
separately for both sides of the boundary. In addition, the eigenvectors are determined 
as a function of k^ and E. 

In the second step, the boundary conditions are applied, finally yielding an energy 
dispersion E{kx)- For hard wall boundary conditions, the wave function is required to 
vanish at the boundary between the HgTe sample and vacuum. For soft wall boundary 
conditions, both the wave function and its first derivative are required to be continuous 
at the boundary between the topologically trivial and the topologically non-trivial HgTe 
samples. In both cases, we require that the solutions are bound to the edge, i.e. that 
the probability densities vanish far away from the boundary (see figure [5]). 
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Figure 5. Normalized probability density at fc^; = as a function of y for soft wall 
boundary conditions separating a topologically non-trivial QW from vacuum (trivial) 
(a) and hard wall boundary conditions separating a topologically non-trivial QW region 
from a topologically non-trivial QW region. 




Figure 6. Energy dispersion with a boundary in x-direction. On the left hand side 
of the boundary there is a topologically non-trivial QW with C'=7.8mcV (black bulk 
energy minimum), while on the right hand side of the boundary there is either vacuum 
(a) or a topologically trivial QW with C""=:6.8meV, indicated by the green bulk energy 
minimum (b). Chiral Majorana edge states obtained numerically as boundary solutions 
to ll]) are shown in red. The red dotted lines are estimations derived from pT|) . The 
parameter values are taken from table [TJ 



Calculation of edge state dispersion for small 

Under certain simplifying assumptions, the edge dispersion of the boundary problem 
can be calculated to first order in k using perturbation theory. We assume that is 
is sufficient to calculate the edge dispersion in the reduced model ([5]) and require that 
quadratic terms oc kl,ky may be neglected (Assumption 1). Consistently, we also match 
only the wave functions and not the derivatives at the boundary. 

We first prove that the boundary conditions are fulfilled for E = 0, k^ = and calculate 
the wave function of the corresponding bound state. 

Imposing the condition E = k^ = in yields four possible values for A = iky-. 

±X±= A ^ • (13) 



'Rl + h^ + 2Roh sin 29 
Depending on the parameter values, the signs of A can be determined as follows under 
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the assumption that \Be\ > |/io| 





5| < /ig + A| 


+A+ > 


+A+ > 


-A_ > 


+A_ > 



The respective eigenvectors are 



|A.) ^ - 



-sign (Be) a/ 1 + fip/B^ 
e^^^sign {BE)^l + fio/BE 



\ 



1 - f^o/Bi 



J 



with e*^ = ( iRo + he 



(14) 



A+) and | — A_) 



/^Rl + h^ + 2RQhsm2e and where |A_), 
can be obtained from |A+) by clockwise translation of the minus sign in front of the 
second entry. 

If B\ > fif + A|; on the left side of the boundary (topologically non-trivial) and 
Be < fJ'r + ^"e '^'^ the right side of the boundary (topologically trivial), the matching 
condition requires that the eigenvectors on the left and right hand sides of the boundary 
are linearly dependent, 

det (lA',), I - A'_), I - a;), I - AD) = (15) 

yielding the solution 



Af 



Q{-y)e 



\i) + Q{y) (cie-^+^l - a;) + cse"^-^! - A^)) j . (16) 

1. The coefficients Ci and C2 



Here A/" is a normalization factor such that J dy\ip{y)\'^ 
are determined by the matching condition. 

This proves that there is a bound state for kr^ = E = 0. The antisymmetry of the 
Hamiltonian ([T]) under particle-hole transformation S^HE = —H, where S = iKayTy, 
has the important consequence that ip{y) for = E = constitutes a Majorana 
excitation, as can be easily confirmed by explicit calculation. 

Now we consider the case of non-zero k^. When kx is small, the fca;-dependent BIA 
and Rashba terms in equation ([5]) can be treated as a perturbation V. The expectation 
value of this perturbation with respect to the kx = E = state gives the energy in 
terms of kx to first order in perturbation theory, 

I A I R^ — ^^^^ 

E = / dyr{y)V^{y) ' ° 



M2 



Br 



(17) 



Af2 



Edge state dispersions 

The resulting edge energy dispersions that were calculated numerically from ([1]) and 
using (frfll are presented and compared to the bulk energy dispersions in figure [6l 

The boundary between topologically trivial and non-trivial regions results in a single 
chiral Majorana bound state of approximately linear energy dispersion that smoothly 
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connects to the occupied and non-occupied bulk bands. The existence of this edge 
state is independent of which boundary condition (soft or hard wall) or which inversion 
asymmetry term (BIA, Rashba or both) is considered. As expected, its existence 
depends on the topological bulk properties on both sides of the boundary - the edge 
state disappears when the Chern numbers on both sides of the boundary are identical. 

The formula f|T7j) describes the numerically obtained energy dispersions of figure [6] 
well, regardless of whether hard wall or soft wall boundary conditions have been 
employed. A small discrepancy between the numerical calculation on the basis of the 
Hamiltonian ([1]) and (ITTI) is mostly due to the limited applicability of Assumption 1. 
Yet, ( ITTll describes the salient features of the edge dispersion well. 

As predicted by f|T71) . the edge state velocity depends most strongly on the 
parameters 9, A, h, Rq, Be, Ae- We have checked that the parameters D, B, C, Ah, Bh 
have comparably little effect on the energy dispersion. 

Interestingly, we find that the group velocity v = of the chiral edge state is 
tunable by the interplay between the Rashba and BIA spin-orbit interactions. The edge 
state is right-moving when the Rashba term dominates and left-moving when the BIA 
term dominates (see figure [7]). This is consistent with the change of the first Chern 
number in (11 II) between Ci = ±1 by means of the index theorem [28]. When applied to 
the hard-wall boundary condition, the index theorem equates the topological charge Ci 
of the bulk with the sign u of the slope of the chiral edge mode. It is therefore easy to 
see that as the Rashba parameter Rq changes from R^ < to Rq > h"^, the first Chern 
number changes from —1 to +1 and the edge state from left to right moving. This could 
be experimentally verified by changing the Rashba parameter, while the intrinsic BIA- 
SOC remains fixed. Another interesting observation is that the group velocity of the 
chiral Majorana edge state depends on the orientation of the edge with respect to the 
crystal axis when both Rashba and BIA SOC terms are present. The absolute values of 
the velocity is smallest along the [110] {6 = n/A) and largest along the [110] {6 = 3tt/4) 
direction. In HgTe-QWs, the edge structures are for experimental reasons prepared in 
parallel to the cleavage planes, and therefore 95 % of the current carrying edges are 
orientated either in [110] or [110] direction [35]. Due to the simultaneous presence of 
Rashba and BIA SOC terms, these two experimentally accessible edge orientations are 
expected to have distinct chiral Majorana edge state properties. 

5. Comparison of chiral Majorana edge states and helical edge states 

As described in the previous section, a HgTe-QW in proximity to a superconductor 
with BIA and a Zeeman field hosts chiral Majorana edge states at its edges when the 
parameters are in the topological nontrivial TSC phase (Ci = ±1). 

On the other hand, it is well-known that a time-reversal invariant HgTe-QW hosts 
helical edge states at its boundary [161 EHl EH] when the mass parameter M is negative 
{B/M > 0). 

In the following, we will compare the chiral Majorana edge states and the helical 
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Figure 7. The edge group velocities function of different parameters. 

Analytic approximations (solid black lines) obtained from ()17p are compared with 
numerical values for hard wall (red triangles) and soft wall boundary conditions (blue 
crosses). This confirms that the edge state only exists when the boundary connects 
topologically trivial and non-trivial regions, compare The Rashba spin-orbit 

coupling is set to zero except in the last line; all other parameters are taken from table[l] 



(Dirac) edge states, assuming that only the BIA-SOC term is present. 

Chiral Majorana edge states and helical edge states at a single boundary are 
mutually exclusive. Firstly, helical edge states are not topologically protected when 
time reversal symmetry is broken by an Zeeman field, which is needed for the emergence 
of chiral Majorana edge states [36]. Secondly, it has recently been shown that even in 
the presence of time reversal symmetry, the introduction of the superconducting term 
Atx adiabatically interconnects the regions that in its absence would differ by their 
topological invariants [37] . 

In the parameter range where there is a chiral Majorana edge state, there are 
therefore no additional helical edge states within the bulk energy gap, see figure [H] (a,b). 
Even if the mass gap M is negative, potential helical edge states are gapped out by the 
proximity of the superconductor. 

For comparison to the chiral Majorana edge states, the helical edge states are shown 
in the absence of both the Zeeman field and the superconductor in figure[8](c,d), where in 
contrast to the standard HgTe the additional BIA term is included. As reported in [38], 
it has a negligible effect, and the velocity of the helical edge states is still well-described 
by [39] 
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for small \kx\. The chiral and helical edge states differ greatly in their velocities, which 
are given by f|T71) and f|T8l) . respectively. The offset is also different: While chiral 
Majorana edge states are confined to the crossing kx = E = hj particle hole symmetry, 
the helical edge states are only restricted to crossing each other somewhere on the E 
axis by time reversal symmetry. 

Helical and chiral edge states have similar widths. While recent calculations for the 
width of helical edge states predict an edge width under 100 nm in a HgTe-QW near 
the Dirac point [39], experimentally confirmed up to a factor of 2 in ^4(1|, the chiral 
Majorana edge states typically extend to about 200nm-300nm in a HgTe-QW for the 
parameters in table [H 




kx[(^im)"''] kx[(^m)"''] 

(a) (b) 




-50 50 -50 50 



kx[(^mr^] kx[(^m)"''] 

(c) (d) 

Figure 8. Energy dispersion with a hard wall boundary in a;-direction for a HgTe- 
QW with topologically nontrivial parameters (a) C = — 7.8meV, M — 8meV 
and (b) C — 7.8 meV, M = — 8meV. This allows chiral Majorana edge states to 
form. All other parameters except the Rashba-SOC, which is set to zero, are taken 
from table [TJ In (c,d), the energy dispersions of a time- reversal invariant HgTe- 
QW in the inverted regime without the presence of the superconductor are shown 
(C = 6.8 meV, M = — 8meV). For better comparison with (b), the degrees of freedom 
have been doubled in (d) by reflection at the Fermi level. 

Summary 

In summary, we have shown that a HgTe/CdTe QW in the proximity to an s-wave 
superconductor and the presence of a Zeeman-splitting can represent a TSC with chiral 
one- dimensional edge states. At kx=E=0, these quasiparticles have the Majorana 
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property of being their own antiparticle. 

Importantly, an external electric field that leads to an additional Rashba-SOC term 
is not required, as long as the intrinsic BIA is large enough. While in most publications, 
BIA is discussed as an unwanted but small perturbation, in this case a high BIA would 
be of interest. Interestingly, if Rashba-SOC and BIA-SOC are both present, the group 
velocity v = \^ of the chiral edge state depends crucially on the interplay between 
both spin-orbit interactions and can be reversed if the Rashba-SOC is tuned from zero 
to a value that exceeds the bulk inversion asymmetry term. We discussed in analytical 
and numerical terms in which parameter regime the TSC exists using the extended four- 
band model of BHZ valid near the topological phase transition between a trivial and a 
non-trivial TI that would exist without the superconducting and Zeeman terms. 

Although many systems hosting chiral Majorana fermions have so far been 
presented [SlIIO], this is the first proposal to switch between one- dimensional Majorana- 
and Dirac edge states in a (HgTe) QW, rather than the surface of a 3D TI, where 
different scenarios to convert between one-dimensional Majorana- and Dirac edge states 
have previously been proposed [HI |l2l |33], and the first time chiral Majorana- and 
helical Dirac edge states appear in the same system. 

The possibility to switch between chiral Majorana- and helical Dirac edge states 
on the same one- dimensional boundary might not only enhance experimental control 
in the detection of Majorana quasiparticles, but could also provide an interesting 
framework that combines dissipationless transport and the non-local storage of quantum 
information. 

Importantly, we would like to point out that our analysis should equally well hold 
for other 2D TIs, like type-II QWs made of InAs/GaSb/AlSb, which have intrinsically 
inverted bandstructures and similar inversion symmetry breaking SOC terms |44j . 
Experimental evidence of the quantum spin Hall effect in these systems has recently 
been given [15], also in connection with superconducting contacts |16]. Further, it has 
been proposed that bilayer HgTe-based QWs in the presence of an interlayer voltage 
allows for an all electrically tunable TI |171 H8] . 
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Appendix 

If the mass gap M is tlie relevant energy scale of tlie problem, leading to only four 
relevant £^l-like energy bands, this allows us to reduce the eight-band Hamiltonian ([1]) 
onto an effective four-band Hamiltonian. 

To this end, we divide the Hamiltonian into three parts, 

H = Ho + Hi + H2 (A.l) 
Ho = {C + Ms,)t, (A.2) 

Hi = Ro ^^'^^\ kcry - kya,^)T, 

+ (A+ + A^s,) T,-{B + D) P + (5+ + a, 
H2 = Ak^s^a^T^ - AkySyT^ + haySyT^. (A.4) 

where for clarity, in contrast to the more general Hamiltonian ([1]), we set ^ = 0. 

In the following, Hq (containing only Sz) is considered as the unperturbed 
Hamiltonian and Hi (containing only Sz) and H2 (containing Sy) as the block-diagonal 
and block-offdiagonal perturbations, respectively. The perturbations are small in the 
sense that the mass gap M is the dominant energy scale, i.e. 

|A±|, \B±\, \h\, \A\\k\, \Ro\\k\, \B\\k\'^, \D\\k\^ < |M| (A.5) 

In order to approximately diagonalize the Hamiltonian H, 

H = e-^He^ (A.6) 

we apply quasi-degenerate perturbation theory ("Lowdin partitioning") |29j . 

The important matrix for the transformation, S", consists of different orders of the 
perturbation parameter, 5* = S^^^ + 0{2). It can be shown that 

H = Ho + Hi + ^[H2, 5(1)] + 0{3) (A.7) 

is block-diagonal up to second order if S^^^ is determined in such a way that 

[Ho,S^'^] = ~H2. (A8) 
For the Hamiltonian ( lA.ip it is straight-forward to sho'ftj^ that 



• ^ .Akx . Aky , . . 

•J = -^^^rrSxOy — i^rrSyOz — i^rrSx (A. 9) 

2M ^ 2M ^ 2M ^ ^ 



H = fi{sz, k)Tz + B{sz)az + A{sz)tx 

hkySzCJyTz + R{sz)kxCTyTz - R{sz)kyaxrz 

where the new parameters 



(A.IO) 



f^{sz, P) = eCky) + (^M(P) + ^ + ^y^ (A.ll) 
under the requirement that the parameters do not depend on x or y. 
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B{s,) = B+ + B_s, (A.12) 
A(s,) = A+ + A_s, (A.13) 

(A.14) 

R(a,) = -Ro^^^^ (A.15) 

have been introduced. The advantage of the transformed Hamiltonian ( lA.lOp is that Sz 
is (approximately) a good quantum number. For further calculations, we neglect the 
S2 = — 1 solutions in the main text ([5]), and retrieve the angular dependence on 9 by 
the simple replacement rules 

hay — )■ h (cos 29ay + sin 26ax) (A. 16) 

ha,^ h (cos 29a.j: - sin 29ay) (A. 17) 

which amount to a rotation in the xy-plane. 
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